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I. INTRODUCTION 


The Wilson line in QCD is extensively used at high energy colliders. In particular, the 
Wilson line is used in the definition of the non-perturbative quantities in QCD such as in 
the definition of the parton distribution function and fragmentation function at high energy 
colliders [1]. There are two types of Wilson lines, 1) light-like Wilson line and 2) non-light- 
like Wilson line. 

In case of light-like Wilson line the parton travels at speed of light and hence it produces 
pure gauge potential at all the time-space points except at the position x transverse to 
the motion of the parton at the time of closest approach 2-4], Since light-like Wilson line 
produces pure gauge potential, it produces zero fields i. e., F^(x) = 0. 

In case of non-light-like Wilson line the parton does not travel at speed of light and hence 
it does not produce pure gauge potential [3-4]. Since non-light-like Wilson line does not 
produce pure gauge potential, it produces non-zero fields i. e., F“ v (x) ^ 0. In addition to 
non-zero F^ u (x), the form of F^(x) can be arbitrary for non-light-like Wilson line. Hence 
one expects fundamental physics differences between light-like Wilson line and non-light-like 
Wilson line. 

One of such physics issue we will study in this paper is the factorization theorem in QCD 
and QED. We will show in this paper, by using path integral formulation of QCD and QED, 
that the factorization theorem is valid for light-like Wilson line but is not valid for non-light- 
like Wilson line. This conclusion is also consistent with Ward identity and Grammer-Yennie 
approximation (see section III for details). 

In case of Grammer-Yennie approximation the polarization sums for both real and virtual 
photons are rearranged into two parts. One of these is called the I\ -polarization sum which 
is proportional to k^k v which corresponds to longitudinal polarization of the photon field. 
This corresponds to pure gauge because the polarization of the pure gauge is proportional 
to its four momentum ki. e., it is longitudinally polarized (which is unphysical) which 
can be gauged away in the sense of factorization [5]. 

In case of Ward identity the gauge transformation part of the photon field is nothing but 
the pure gauge. Since the polarization of the pure gauge is longitudinally polarized (which 
is unphysical) it can be gauged away in the sense of factorization (gauge transformation). 

Since the light-like Wilson line produces pure gauge and the non-light-like Wilson line 
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does not produce pure gauge, one finds that the Ward identity and Grammer-Yennie ap¬ 
proximation is consistent with the fact that the factorization theorem is valid for light-like 
Wilson line but is not valid for non-light-like Wilson line. 

Note that for non-light-like Wilson line there is no collincar divergence in the Eikonal 
approximation [see eq. flD) below], ffence for non-light-like Wilson line one does not have 
to worry about collinear divergences. In case of non-light-like Wilson line there can be soft 
divergences which are not factorized which we will prove in this paper. 

Note that, in some of the studies, the non-light-like Wilson line is used in the definition 
of the transverse momentum dependent parton distribution function (TMD PDF) at high 
energy colliders jl, 6]. These studies involve diagrammatic calculation at one loop level using 
perturbative QCD. However, the parton distribution function and fragmentation function 
are non-perturbative quantities in QCD. Hence if one uses the perturbation theory to study 
their properties l|, 6], then one may end up finding wrong results. In general, the non- 
perturbative phenomena may be impossible to understand by perturbation theory, regardless 
of how many orders of perturbation theory one uses. Take for example, the non-perturbative 
function 


/0) = e" 


n = positive integer. 


( 1 ) 


The Taylor series at x = 0 for this function f(x) is exactly zero to all orders in perturbation 
theory, but the function is non-zero if x ^ 0. 

Hence the diagrammatic method using perturbative QCD may not be always sufficient 
to prove factorization of soft-collinear divergences of parton distribution function and frag¬ 
mentation function at high energy colliders, regardless of how many orders of perturbation 
theory one uses, because parton distribution function and fragmentation function are non- 
perturbative quantities in QCD. 

On the other hand the path integral method of QCD can be used to study non- 
perturbative QCD. The path integral method of QCD is valid at all order in coupling 
constant. Hence we find that the path integral method of QCD is very useful to prove 
factorization of soft-collinear divergences of non-perturbative quantities in QCD such as the 
parton distribution function and fragmentation function at all order in coupling constant at 
high energy colliders. 

In this paper we prove, by using path integral method of QCD, that the factorization 
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theorem in QCD and QED is valid for light-like Wilson line but is not valid for non-light- 


like Wilson 


studies jl, 6 1 which used the non-light-like Wilson line in the definition of the (transverse 
momentum dependent) parton distribution function and fragmentation function at high 
energy colliders. This is in conformation with the finding in [7, 8 ] which proved factorization 


ine. This implies that the factorization theorem is violated in all the previous 


0 


and 


theorem in NRQCD heavy quarkonium production in case of light-like Wilson line 
the violation of factorization theorem in NRQCD heavy quarkonium production in case of 
non-light-like Wilson line 1 8 ]. In case of massive Wilson line in QCD the color transfer 

n 

occurs and the factorization breaks down [ 8 J. The simple physics reason behind this is the 
following. The light-like Wilson line produces pure gauge potential which gives F“ v (x) = 0 
which can be gauged away in the sense of factorization because pure gauge corresponds to 
unphysical longitudinal polarization and hence factorization theorem works (see sections IV 
and VII for details). The non-light-like Wilson line does not produce pure gauge potential 
{F“ v (x) 7 ^ 0), the form of which can be arbitrary which can not be gauged away in the sense 
of factorization because non-pure gauge does not correspond to unphysical longitudinal 
polarization and hence the non-light-like Wilson line spoils the factorization (see sections V 
and VIII for details). 

The paper is organized as follows. In section II we briefly discuss soft-collinear divergence 
and light-like Wilson line in QED. In section III we discuss Ward Identity, Grammer-Yennie 
approximation, longitudinal polarization and light-like Wilson line. In section IV we discuss 
the proof of factorization of soft-collinear divergence in QED by using light-like Wilson line. 
In section V we prove that the soft-collinear divergence in QED is not factorized when non- 
light-like Wilson line is used. In section VI we discuss soft-collinear divergence and light-like 
Wilson line in QCD. In section VII we discuss the proof of factorization of soft-collinear 
divergence in QCD by using light-like Wilson line. In section VIII we prove that the soft- 
collinear divergence in QCD is not factorized when non-light-like Wilson line is used. Section 
IX contains conclusions. 
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II. SOFT-COLLINEAR DIVERGENCES AND LIGHT-LIKE WILSON LINE IN 
QED 


As mentioned earlier in order to study factorization of soft-collinear divergences the Wil¬ 
son line is used in the definition of parton distribution function and fragmentation function 
which also make them gauge invariant. Before proceeding to the issue of gauge invariance 
and the Wilson line for soft-collinear divergences in QCD let us first discuss the correspond¬ 
ing situation in QED. The gauge transformation of the Dirac held of the electron in QED 
is given by 

if/(x) = e ieu[x) ^{x). (2) 


Hence we can expect to address the issue of gauge invariance and factorization of soft- 
collinear divergences in QED simultaneously if we can relate the c o(x) to the photon held 
A^(x). This is easily done by using Eikonal Feynman rules in QED for soft-collinear diver¬ 
gences when light-like Wilson line is used. This can be seen as follows. 

The Eikonal propagator times the Eikonal vertex for a photon with four momentum k^ 
interacting with an electron moving with four momentum p 1 ' in the limit k «< p is given 


by 


QS, 0E- 


16] 


e — TA~~- 3 

p ■ k + le 

In QED the soft divergence arises only from the emission of a photon for which all com¬ 
ponents of the four-momentum k^ are small (k —> 0) which is evident from eq. ([3]). From 
eq. ([3]) one also finds that when 0 < k <<< p and p is parallel to k one may find collinear 
divergence. 

However, since 


_> / 777 , 

p . k = p 0 k 0 - p- k = \p\\k\(Jl + —— cos#) 


( 4 ) 


one finds that the collinear divergence does not appear in QED because of the non-vanishing 
mass of the electron, i. e., m ^ 0. From eqs. (j3J) and (HJ) one finds that the collinear 
divergence appears only when m — 0 and 9 = 0 where 9 is the angle between p and k. 
Since gluons are massless and the massless gluons interact with each other one Ends that 
the collinear divergence appears in QCD. Since a massless particle is always light-like one 
finds that the soft-collinear divergences can be described by light-like Wilson line in QCD. 
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For light-like electron we find from eq. fl3]) 


e —--- = e 

p ■ k + ie 


l > 1 

l ■ k + ie 


(5) 


where is the light-like four-velocity (|/j=l) of the electron. Note that when we say the 
” light-like electron” we mean the electron that is traveling at its highest speed which is 
arbitrarily close to the speed of light (|/| ~ 1 ) as it can not travel exactly at speed of light 
because it has finite mass even if the mass of the electron is very small. Hence we find 
that if is light-like four velocity then the soft-collincar divergence can be described by the 
Eikonal Feynman rule as given by eq. (J5]). 

From eq. (E]) we find 


d 4 k l ■ A(k) 

(27r) 4 l ■ k + ie 


roo r rj^k r°° 

/ dX / — —e il - kX l • A[k) = ie. / dXl ■ A(IX) 
Jo J (27T ) 4 JO 


= — ei 


(2t r 


( 6 ) 


where the photon held A^[x) and its Fourier transform A^(k) are related by 

r d 4 k 
A “<*> = / 


A k-x 


(7) 


Now consider the corresponding Feynman diagram for the soft-collinear divergences in QED 
due to exchange of two soft-collinear photons of four-momenta k± and k%. The corresponding 
Eikonal contribution due to two soft-collinear photons exchange is analogously given by 


2 r d% d 4 k 2 l • A(k 2 )l • A(k i) 

J (27 r) 4 ( 27 t ) 4 (/ • (ki + k 2 ) + ie)(l ■ ki + ie) 

r oo poo 

= e 2 i 2 / dX 2 / dXA ■ A{lX 2 )l ■ A{1 Ai) 

Jo J a 2 

/*oo poo 

= — r dX 2 dXA ■ A(lX 2 )l ■ ^4(ZAi). (8) 

2! Jo Jo 


Extending this calculation up to infinite number of soft-collinear photons we find that the 
Eikonal contribution for the soft-collincar divergences due to soft-collinear photons exchange 
with the light-like electron in QED is given by the exponential 


e ie ir dXl ' A dA 


(9) 


where is the light-like four velocity of the electron. The Wilson line in QED is given by 


ie dx^Auix) 
e Jx i 


( 10 ) 
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When H M (x) = A* 1 (l\) as in eq. (J9]) then one finds from eq. (flUl) that the light-like Wilson 
line in QED for soft-collinear divergences is given by [hzj 


Je£f dx^A^x) _ —ie /“ d\l-A{x f +l\) ie d\l-A( Xi +lX) 


( 11 ) 


Note that a light-like electron traveling with light-like fonr-velocity F produces U(l) pure 
gauge potential A^(x) at all the time-space position x' 1 except at the position x perpendicular 
to the direction of motion of the electron (/ • x — 0) at the time of closest approach ffl. 
When A 1 ' 1 (x) = A^(\l) as in eq. © we find l ■ x — XI ■ l — A ^ 0 which implies that the 
light-like Wilson line finds the photon held A^(x) in eq. (EJ) as the U(l) pure gauge. The 
U(l) pure gauge is given by 

A fl (x ) = d tl co(x) ( 12 ) 

which gives from eq. (TTT1) the light-like Wilson line in QED for soft-collinear divergences 


e ieu(x f ) e -ieu(xt) = £ ie £. f = ^ d \l-A(x f +l\) £ ie d\l-A( Xi +l\) 


(13) 


which depends only on end points xf and x j but is independent of the path. The path 
independence can also be found from Stokes theorem because for pure gauge 


F^(x) = d^A v (x) - d v A»(x) = 0 


which gives from Stokes theorem 


o ie j> c dx/ 1 A /J ,(x) _ ie J s dx^dx u (x) ^ 


(14) 


(15) 


where C is a closed path and S is the surface enclosing C. Now considering two different 
paths L and M with common end points xf and x^ we find 

e ie § c dx ^ A A x ) — e ie f L dx ll Afj,(x)—ie f M dx»A ll (x) _ ^ (qg\ 


which implies that 


ie J x f dx * 1 A^x) 


(17) 


depends only on end points xf and x^ but is independent of path which can also be seen 
from eq. (TT3|) . Hence from eq. (TTH]1 we find that the abelian phase or the gauge link in QED 
is given by 

e ~ie f 0 °° d\l-A(x+l\) _ £ ieuj(x) Qg\ 
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Eq. (jT8|) can also be directly verified as follows. When A tl {x) = d^uj^x) is U(l) pure gauge 
as given by eq. (fl2|) we find by using contour integration that 


e -ie f 0 °° d\l-A(x+l\) = e -iedXe™ l-A(x) = g -ie f 0 °° d\e™ a f ^e^l-A(k) 

= p - ie f<0zf O °° dXeiX ‘ k = e^l-A(k) = ie _l_ t . A(x) = ietJj(x) 


(19) 


which reproduces eq. m 

From eqs. (121) and (TT8h one expects that the gauge invariance and factorization of soft- 
collinear divergences in QED can be explained simultaneously if the light-like Wilson line is 
used. 


III. WARD IDENTITY, GRAMMER-YENNIE APPROXIMATION, LONGITU¬ 
DINAL POLARIZATION AND LIGHT-LIKE WILSON LINE 


Note that, as mentioned above, in classical electrodynamics the assertion 


ihat the gauge 


field that is produced by a highly relativistic particle is a pure gauge ( 2 - 4]. I 11 order to 
study factorization of infrared divergences by using the background field method of QED, 
the soft photon cloud traversed by the electron is represented by the background field 
( 5 )] where one represents the (hard) quantum photon field by Q^(x). 

In order to see how, in the quantum field theory (QED), the pure gauge property holds 
for the case in which the photon momentum is soft can be seen from the Grammer-Yennie 
approximation [lfj, see also jfj. I 11 the Grammer-Yennie approximation, if the photon is 
emitted from an incoming line, then for soft-divergence one replaces the polarization of the 
gauge field [see eq. (3.2) of jl8[] by 


e“(k) = PW - vQH] + U ' 


l ■ k 


l-k 


( 20 ) 


which, when used in the Feynman rule for the eikonal line [that is given in eq. 
in the equation before eq. (3.2) of Q], gives 


Z-e(fc) rl-e(k) 


l-k 


l ■ k 


l-e(k) ] | l-e(k) Q | l-e(k) 


l-k 


l ■ k 


l-k 


above or 


( 21 ) 


Hence one finds that the first term [e M ( k ) — k in the right hand side of eq. (I2(7|i can not 
lead to soft divergence. The soft divergence comes from the last term k^ 1 -^- in the right 
hand side of eq. (l20|h Hence from the Grammer-Yennie approximation one finds that, as 
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far as soft divergence is concerned, the polarization of the gauge held is proportional to k^, 
i. e., it is longitudinally polarized. Since the polarization of the pure gauge is proportional 
to k one finds that the pure gauge property holds for the case in which the gauge held 
momentum is soft. Since longitudinal polarizations are unphysical they can be gauged 
away, in the sense of factorization fj]. Hence one hnds that the factorization theorem for 
soft/infrared divergences can be studied by using pure gauge. 

The pure gauge property is also evident from the Ward-Identity which can be seen when 
one replaces e M (/c) in the QED amplitude 

M = e»{k)M^k) (22) 

by 

e >J '(k) —> k M Pure gauge (23) 

then one gets 


k^M^k) = 0 (24) 

which eliminates longitudinal polarization of the photon in the physical process. This is the 
statement of gauge transformation which is equivalent to saying that the gauge transforma¬ 
tion part is noting but the pure gauge. Hence the statement that unphysical longitudinal 
polarization can be gauged away by using pure gauge, in the sense of factorization B , is 
consistent with Ward-identity. 

The above analysis in the quantum held theory is also valid for collinear divergence in 
the Eikonal approximation for light-like Wilson line case because in the collinear limit the 
momentum k 11 of the gauge held is proportional to the four velocity of the light-like Wilson 
line Since the polarization of the pure gauge is proportional to one hnds that the 
pure gauge property holds for the case in which the gauge held momentum is collinear to 
the momentum of the light-like Wilson line. 

Note that for non-light-like Wilson line there is no collinear divergence in the Eikonal 
approximation which can be seen from eq. (|4j). Hence for non-light-like Wilson line one does 
not have to worry about collinear divergences. 

In case of non-light-like Wilson line there can be soft divergences which are not factorized 
which we will prove in this paper. This is consistent with the fact that non-light-like Wilson 
line can not produce pure gauge j^-4|. 
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Hence, since the light-like Wilson line produces pure gauge and the non-light-likc Wilson 
line does not produce pure gauge [2-4], one finds that the Ward identity and Grammer- 
Yennie approximation is consistent with the fact that the factorization theorem is valid for 
light-like Wilson line but is not valid for non-light-like Wilson line. 


IV. PROOF OF FACTORIZATION OF SOFT-COLLINEAR DIVERGENCES IN 
QED WITH LIGHT-LIKE WILSON LINE 


Note that the gauge invariant greens function in QED 

rx 2 

G(xi, x 2 ) =< ^(^ 2 ) x exp[ie / dx^A^x)] x 1 ) > 


(25) 


/Xi 


in the presence of background field A IJ, (x) was formulated by Schwinger long time ago 19]. 
When this background field A^(x) is replaced by the U(l) pure gauge background field as 
given by eq. (1T2|) [which corresponds to the light-like Wilson line in eq. f[I5|) ] then the 
path integral method of QED can be used to prove gauge invariance and factorization of 
soft-collinear divergences in QED simultaneously. 


A. Proof of Factorization Theorem in QED With Light-Like Wilson Line in Co¬ 
variant Gauge 

The generating functional in QED in the path integral formulation in covariant gauge is 
given by 

Z[J, 77 , fj\ = j[dQ}[d$}[dil>\ e i J d4x[ -^ [Q] -^ d ^ )2+ ^ i ^-™+^QM+J-Q+v^+M ( 2 6) 

where J M (x) is the external source for the quantum photon field Q^(x) and fj(x) is the 
external source for the Dirac field ip{x) of the electron and 

F“[Q] = di‘qr(x)-ff'cr(x), FlM = (27) 

The correlation function of the type < O|'0(cc 2 )'0(2a)|O > in QED in covariant gauge is given 
by 0 

< 0|tK.t 2 )^Od)|0 >= J[dQ][d$][di/;] 

(28) 


10 






where suppression of the factor is understood as it will cancel in the final result (see eq. 

W 

The generating functional in the background field method of QED in covariant gauge is 
given by 

Z[A,J,n,ij] = J{dQ ][#][#] 

(29) 

where we have used the notation Q^(x) for the quantum photon field and A^(x) for the 
background field. The correlation function of the type < 0|'?/>(.T2)t/ ; (^i)|0 > in the background 
field method of QED in covariant gauge is given by jfj 

< 0|Vi(x 2 )i/’(ii)|0 >.4= J[dQ] [<$][#] 

(30) 

where suppression of the factor is understood as it will cancel in the final result (see eq. 

(I39]) h We write eq. (130]) as 

< 0|qx 2 )</>Ui)|0 >/i= f[dQ][di>'][dili'} ^'( X2 )^'( Xl )e i /' <4l [-J F 2»l«l-i:( fl i-0') 1 +»'[n , ‘9»-”+«7'‘(/i+OW*'] 

(31) 


because the change of integration variable from unprimed to primed variable does not change 
the value of the integration. 

The light-like Wilson line in the background field method of QED is given by, see eqs. 
(HU) and HE]) 

e ie fo dx'i*A^x') ( 32 ) 

where A^{x) is the U(l) pure gauge potential as given by eq. (fT2j) . Hence for light-like 
Wilson line we find 

[ dx' v A v (x') = / dr — ^ ^ — = utix) — a;(0), xYl) = x, x'(0 ) = 0 (33) 

Jo Jo dr 

which gives 

<9 m [ f dx' u A v {x')] = d^ui{x) = A^x). (34) 

J 0 
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From eq. (TT8|) let us define 

= e~ ie fo°° dXlA ( x+lx ^( x ) = e <ew (*V(®)- (35) 

From eq. (TTU]) we find for light-like Wilson line [H M (x) = d^uj^x)] that 

d\l ■ A(x + IX)] = ~d,,{j—Q 1 ■ A ( x )] = -d fl u(x) = -A^x). (36) 

Hence when the background held A ll {x) is the U(l) pure gauge as given by eq. (1T2T) . 
which is the case for light-like Wilson line, we find from eqs. (1351) . (1361) and (T311) that 

< 0|^(X 2 )^(X!)|0 > A = J[dQ] [#][#] ^(^ 2 )e ie -/o°° rfAZ.^( a , 2 +iA) e -ie/ 0 °° 


Xe* / d4x l-j F ^AQ]-^;(d )1 ,Q ll ) 2 +4>[i'y^d l _ l -m.+e'y> I -Q IJ ,]i'} (gy) 

which gives by using eq. (USD 

< 0\i>(x 2 )e e ^ x i dx ‘ 4ai(x V(^i)|0 >a= I [dQ] [dip] [dip] i>(x 2 )i>(x 1 ) 

Hence from eqs. (138]) and (128|) we find that 

< O|'0(x2)'0(^i)|O >=< 0|^(x 2 )e* e -^i dx A ^ x \{x i)|0 (39) 

where < O|'0(x 2 )e* e -^ : i dx A ^ xp{xi) |0 is gauge invariant. By using eq. (1T3T) we find from 
eq. (1391) that 

< ^(x 2 )^(xi) >=< ${x 2 ) e~ ie fo° dXl ' A ( X2 + lx ) e ie ir dXl ' A ^ xl+lx ^ ^{xi) > A (40) 


which proves factorization of soft-collincar divergences in QED at all order in coupling 
constant in covariant gauge when light-like Wilson line is used. The eq. (1401) was first 
proved by Tucci in jfj]. We have included its derivation in this section because we will follow 
the similar path integral technique in the next section to prove that the factorization theorem 
is violated in QED if we replace the light-like Wilson line by the non-light-like Wilson line. 

B. Proof of Factorization Theorem in QED With Light-Like Wilson Line in Gen¬ 
eral Axial Gauge 

The factorization of soft-collinear divergences in QED at all order in coupling constant 
in covariant gauge using light-like Wilson line which we proved in the previous sub-section 
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can also be proved in other gauges. The generating functional in QED in the path integral 
formulation in general axial gauge is given by [ 2 O, 21] 


z[J,v,v] = 


e i f d 4 x[-jF^[Q]-^(r] M Q' J -) 2 +ipli'y' i d M -m+e'y> J -Q ti ]ip+J-Q+rri/)+i/>ri] 


where r is an arbitrary but constant four vector 


V 2 < 0, 

axial gauge 

7? 2 = 0, 

light — cone gauge 

rf > 0, 

temporal gauge. 


(42) 

The correlation function of the type < O|'i/’(;r 2 )'0(aq)|O > in QED in general axial gauge is 
given by 


< O|'0(£ 2 )V , Qh)|O >= 


-0(x 2 )^(x i)e i / d4a; [ _ i- F ^[ < 31“^DM < 9 M ) 2 +V’[n' i ^- m +e7 

(43) 


The generating functional in the background field method of QED in general axial gauge 
is given by 


Z[A,J,V,V] = 


0 i f d A x[-~F%„[Q]--±(ii l _ l Qi 1 ) 2 +^[i'yVd fJ ,-m+e'Y ll (A+Q) ll .]ip+J-Q+r)i’+ipri\ 


(44) 


The correlation function of the type < O|'0(a; 2 )^(a;i)|O > in the background field method of 
QED in general axial gauge is given by 


< O|V ; (^2)'0(^l) |0 >A= 

We write eq. (1451) as 
< 0|^(.t 2 )^(^i)|0 >a= 


'll>(x 2 )'tp(xi)e i / dix \r\ F »AQ}-^(v»Q*) 2 +^[i^d»- r n+e')i J '( A +Q)v 


(45) 




( 46 ) 


because the change of integration variable from unprimed to primed variable does not change 
the value of the integration. 
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Hence when the background field A^{;x) is the U(l) pure gauge as given by eq. (jT 2 T) . 
which is the case for light-like Wilson line, we find from eqs. ([35]) . (I36j) and (146]) that 


< O|'0(x 2 )'0(xi)|O >a— j [dQ][#][#] ip{x2)e ie S™ d ^ A ( x ^) e ~ ie ir dXl - A ^ +lx ^( Xl ) 
xe i f di A~j F j}AQ]-^(vvQ ,, ') 2 +'$[n ,I 'diJ.-™+e'y fi QM (470 

which gives by using eq. (HD 

< 0\4>(x 2 )e e ^i dxliA ^ilj(x i )|0 >^= J [dQ] [df] [df\ x 2 )f>(xi ) 

Xe > / d*x{-\F*„[Q]-±^Q») 2 +$[i^d IJ ,-rn+e'i»Q IJ \'il)\' 

Hence from eqs. (148D and (|43l) we find that 

0\^(x 2 )'ip(xi)\0 >=< 0\ip(x 2 )e e ^ dx Atl ^^(x i)|0 >a 


< 


(49) 


where < 0\ip(x 2 )e e ^ dx A ^ ip{x i)1 0 >a is gauge invariant. By using eq. (TT3|) we find from 


eq. 


that 


< $(x 2 )f’(x i) > = < $(x 2 ) e -ief 0 °°dMA(x2+l* ^ief" dM-Afa+lXty^ >a 


(50) 


which proves factorization of soft-collinear divergences in QED at all order in coupling 
constant in general axial gauge when light-like Wilson line is used. 

C. Proof of Factorization Theorem in QED With Light-Like Wilson Line in Light- 
Cone Gauge 


The light-cone gauge corresponds to 20 22 


V ■ Q = 0, 


rf = 0 


(51) 


which is already covered by eqs. (HI]) and (1421) where the corresponding gauge fixing term 
is given by -i( rj^Q M ) 2 . In the light-cone coordinate system the light-cone gauge 22 ] 


Q + — 0 


(52) 


corresponds to 


rf = (V + ,V ,V±) = ( 0 , 1 , 0 ) 


(53) 
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which covers rj ■ Q = 0 and rj 2 = 0 situation in eq. ()5T]) . 

Since eq. (I50j) is valid for general axial gauges [as given by eq. (jj2j) ] it is also valid for 
the light cone gauge rj 2 = 0. Hence from eq. (HDD we find in the light cone gauge that 

< feMxr) >=< e~ ie -fo° d\i-A( X2 +i\) e ie f™ dXl ' A ^ xl+lx '>^(^ Xl ) > A (54) 

which proves factorization of soft-collincar divergences in QED at all order in coupling 
constant in light cone gauge when light-like Wilson line is used. 


D. Proof of Factorization Theorem in QED With Light-Like Wilson Line in Gen¬ 
eral Non-Covariant Gauge 


The generating functional in QED in the path integral formulation in general non¬ 
covariant gauge is given by 23], 24] 

Z[J, v ,rj\ = f (dQ) [<#][*/>] e ‘f A <-i F 2~l Q 


(55) 


where is an arbitrary but constant four vector. The correlation function of the type 
< 0\'ip(x2)'il>{xi) |0 > in QED in general non-covariant gauge is given by 


< 0\i>(x 2 )'ip{x 1 )\0 >= / [dQ][d4>][d4’\ i>{x 2 )i>(x i)e 


(56) 


The generating functional in the background field method of QED in general non¬ 
covariant gauge is given by 

Z[A, J, i,, rj\ = /[dQ](<#][*/’] e ‘S d4 ’'-FSM-i^^Q,r+iiirs,-m+er(A+Q)MWQ+iM.+H^ 


(57) 


The correlation function of the type < 0|^(x2)V ; (^i)|0 > in the background field method of 
QED in general non-covariant gauge is given by 


< 0\'i/j(x 2 )i>(x 1 )\0 > A = J[dQ][drp][di/}\ i>(x 2 )^(x 1 ) 

e * / d4 x[-\F^[Q]-^(^-d fl Q^ 2 +'lp[i'y^a fl -m+e'y> J -(A+Q) l j\il)\ 


(58) 


We write eq. (j58l) as 

< O|^(x 2 )'0(^i)|O >a= J[dQ][d^i'][dij}'] ^ , (x 2 )'ip'(x 1 ) 

f d 4 x[-\F^[Q]-^(^-dij.Q 1 ,) 2 +4>'[i'y ll d l j,-m+e^ IJ -(A+Q) IJ ,]ip'] 


( 59 ) 
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because the change of integration variable from unprimed to primed variable does not change 
the value of the integration. 

Hence when the background field A^(x) is the U(l) pure gauge as given by eq. flT2j) . 
which is the case for light-like Wilson line, we find from eqs. (I35[) . (I36p and (1551) that 


< 0|^(x 2 )^(a;i)|0 > A = J[dQ][d#iP] i){x 2 )e ie ^ dxi-A( X2 +ix) e -ie f 0 °° dXl - A ^ 1+lx ^(x 1 ) 


xe 


i f d 4 x[-^F^ l/ [Q]-^(^-d ll Q l/ ) 2 +ip[i'y ll d l j,-m+e'y' 1 Q l j,]'il;} 


(60) 


which gives by using eq. (1131) 

< 0\i’(x 2 )e e ^i >A= J [<i-0][c?-0] ${x2)^{xi) 

i f d 4 x[-\F 2 „[Q]--^C^-d^Q^) 2 +‘ijj[i'y' i d^-m+e'yl J -Q^]ip\ 


xe 


Hence from eqs. (15D and (156]) we find that 

< 0\’4i(x2)i^(xj)\0 >=< 0|^(.T 2 )e ie ^i dxl ‘ A ’‘ i ^ijj(xi )|0 > A 


(61) 


(62) 


— ie f X2 dx^A (x) 

where < 0 \'ip(x 2 )e Jx i M ^(xi )|0 is gauge invariant. By using eq. (1T31) we find from 
eq. ( 162 |) that 

n oo p OO 

< ^(x 2 )tP( X i ) > = < ${X 2 ) e~ ie JO d ^ A (x 2 +lX) e ieJ 0 d\l-A( X1 +lX)^ x j >A ( 63 ) 

which proves factorization of soft-collinear divergences in QED at all order in coupling 
constant in general non-covariant gauge when light-like Wilson line is used. 

E. Proof of Factorization Theorem in QED With Light-Like Wilson Line in Gen¬ 
eral Coulomb Gauge 

The generating functional in QED in the path integral formulation in general Coulomb 
gauge is given by 23i. 1241 


Z[J, 77, 77 ] = j [dQ] [ch/j\ [d'ljj] e 4 f d4 A-i F ZAQ]-^([9^- !2 ^]dvQA 2 + 4 P[i^di*-™+e^Q»]4>+J-Q+v^+i’v] 


(64) 


where 


71 " = (1,0, 0,0). 


(65) 
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The correlation function of the type < O|'0(x2)'0(ah)|O > in QED in general Coulomb gauge 
is given by 


< 0\'i/j(x 2 )i>(x 1 )\0 >= j[dQ}[di>][d^\ 

j J ^x[-^[Q]-^([g^-^]d^) 2 +^[ird^-m+erQM 


( 66 ) 


The generating functional in the background field method of QED in general Coulomb 
gauge is given by 

Z[A, J. rj.fj} = j [dQ\ [dip] [dijj] e */ d4 |r [-l i; ’^[ ( 3]-^(b MI/ - 2 y5-] 9 M < 9D 2 +V i [*7 M 9M- m + e 7' 1 (^+Q)d^+J''Q+»7V'+'0»7]_ 


(67) 


The correlation function of the type < Q\'^>{x 2 )'il:{xi)\Q > in the background field method of 
QED in general Coulomb gauge is given by 


< 0| 'il;(x 2 )ip(x 1 )\0 > A = J[dQ][d^][d^\ ip(x 2 )'ifj(x 1 ) 
i J d 4 x[-^F^[Q]-^([gi J ' , '-^^]d IJ ,Q v ) 2 +ip[i'yi 1 -d l ,-m+e'yi J -(A+Q) IJ 


( 68 ) 


Hence by replacing ^- > [g^ — everywhere in the derivations in the previous sub¬ 

section we find 


P oo p oo 

< ^{x 2 )^{ Xl ) > = < ${x 2 ) e~ ie Jo dM.A( X 2 +l\) e ieJ 0 dXl-A(x 1 + l\)^ x j >a ( 69 ) 


which proves factorization of soft-collincar divergences in QED at all order in coupling 
constant in general Coulomb gauge when light-like Wilson line is used. 


V. VIOLATION OF FACTORIZATION THEOREM OF SOFT-COLLINEAR DI¬ 
VERGENCES IN QED WITH NON-LIGHT-LIKE WILSON LINE 

The main reason why factorization theorem of soft-collinear divergences in QED is vi¬ 
olated when the non-light-like Wilson line is used can be easily seen from eq. (]34l) which 
is only valid for the light-like Wilson line but not valid for the non-light-like Wilson line. 
When the Wilson line is non-light-like then in addition to A^x) [like in the right hand side 
of eq. ([31]) ] one will also have a term containing non-zero F^^x), the form of which can be 
arbitrary, which spoils the factorization theorem. 
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For the non-light-like Wilson line we find 


C x 3 /*i dr ,u (r) 

J ^ dx ,v A v {x')] = -^[l dr —— -A v {x\t))], x'(0) = 0, x'(l) = x (70) 


which gives 


rl ,<([£3111 


=/;<M®w.' W )+ a i:^ ))aT a l r) ^ T) ]. (7i) 


dr 


dx' 5 {r) dx /l dr 


Integrating by parts we find 


< ^ - jf 


rl j, dA v {x'{r )) dx ,s (r) dx ,u (r) 
3 dx ,5 (r ) dr 


(72) 


which gives 


da/M^a/)] = A^rr) drF u5 (x\T))^ d -^l. (73) 

Note that for pure gauge we have 


7v((c) = 0 


(74) 


which when used in eq. in (173]) reproduces eq. (134|) which is valid for the light-like Wilson 
line. Hence one finds that for the non-light-like Wilson line, it is the term containing non¬ 
zero F /W (x) in the right hand side of eq. (173]) . the form of which can be arbitrary, which 
spoils the factorization theorem of soft-collinear divergences in QED. 


A. Violation of Factorization Theorem in QED With Non-Light-Like Wilson Line 
in Covariant Gauge 


Note that eq. (130]) is valid for any arbitrary background held A^(x). Hence we start from 
eq. (130]) by following exactly the same procedure that was followed for light-like Wilson line 
case. By changing the integration variable from unprimed variable to primed variable we 
find from eq. (130]) 


< 0\^(x 2 )i>(x 1 )\0 >A= 


J[dQ] [dip'][dip'] 


i F ZAQ]-^;( d ixQ >J ') 2 +4’'[ i 'y' 1 d )J ,-m+e^(A+Q) IJ ,]‘ip'] 


(75) 
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This is because a change of integration variable from unprimed to primed variable does 
not change the value of the integration. For the non-light-like Wilson line let us define the 
primed variable in analogous to the light-like Wilson line case [similar to eq. (1351) ] 

ii/{x) = e- ie ^ dXv ^ x+vX) ^{x) (76) 


where tU is the non-light-like four velocity and A^(x) is not the U(l) pure gauge. Using 
contour integration we find 

d 4 k 


dXv ■ A(x + vX) = I dXe Xv d v ■ A{x) = 


dXe 


Xvd 


e ik ' x v • A(k) 


(2tt) 


d 4 k r c 
(27t) 4 Jo 


dXe 


iXv-k Ak-x 


e lk ' x v ■ A(k) = - 


d A k 


(27r) 4 iv ■ k 


1 e ik ' x v • A(k) = -^ v ■ A(x) 


v ■ d 


(77) 


which gives [since A^(x) ^ d ll uj(x)\ 

roo 2 

J o dXv ■ A(x + vX) = -d/j .[—q v ■ A{x)\ ^ -A^x). (78) 

Note that when A^(x) is the U(l) pure gauge [A ll [x) = d^u ;(x)] we find from eq. (1361) that 
<9 /t [/ 0 °° dXv ■ A{x + uA)] = —A^(x) which precisely cancels the A 11 present in the lagrangian 
density in the background field method of QED [see the derivation of eq. (1371) from eq. 
(13T1) ]. This is the reason why the soft-collinear divergences are factorized in QED when 
A^(x) is the U(l) pure gauge. However, when A lx (x) is not the U(l) pure gauge then such a 
cancelation of the AL M (x) in the lagrangian density in the background field method of QED 
does not happen [see eq. (178]) ] which is why the soft-collinear divergences are not factorized 
in QED when H M (x) is not the U(l) pure gauge. 

From eq. (1781) we find for the non-light-like Wilson line [A^A) ^ d^ l u(x)\ that 

e ~ie d\v-A(x+v\) e ieu{x) (yg) 


where A IJ, (x) is not the U(l) pure gauge. 

Hence we hnd that the equality relation holds for light-like Wilson line, as in eq. (jT8|) 
where A fl (x) is the U(l) pure gauge, whereas the equality relation does not hold for non- 
light-like Wilson line, as in eq. (1791) where A fl (x) is not the U(l) pure gauge. By using eqs. 
(1791) and (1771) in (1751) we hnd 

< 0|^(x 2 )^(.xi)|0 >A= f [dQ][d$][dl(>] ^(072)e ie /o°° rfA^-^(a7 2 +,7A) e -*e/J 30 +^A)^(^ 1 ) 
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which gives 


< 0\ip(x2)e~ ie f 3°° dXv ' A ( X2 + uA ) e ie fo° dXv ' A ( x i + vX ) -0( Xl )1 0 >^ 4 = J [dQ][dif>][dil>] ^(x 2 )^(a:i) 

X g® f d4x [-i F 2v[Q]-j^(dvQ IJ ') 2 +$[ i 7' J ' d v- m + e 7 ,J '( A +Q)v- e 7 ,1 d /1 [^Q v-A(x)}]ili\. (gl) 

Hence from eqs. (1781) . (j 8 ll) and (j28l) we find 

< O| , 0 (x 2 )' 0 (xi)|O >t^< 01 ^(x 2 ) e~ ie fo° dXv ' A ( X2 + uA ) e ie / 0 °° ^-^Oi +t,A )^( Xl ) 10 > A (82) 

which implies the violation of the factorization theorem for soft-collinear divergences in 
QED at all order in coupling constant when the non-light-like Wilson line is used in covariant 
gauge. Note that the equality relation holds for light-like Wilson line, as in eq. (HOll . proving 
the factorization theorem of soft-collinear divergences in QED when light-like Wilson line is 
used. 

Hence by using path integral formulation of QED we have proved that the factorization 
theorem is valid for light-like Wilson line [see eq. (I40j) ] but the factorization theorem is not 
valid for non-light-like Wilson line [see eq. (182]) ]. 

Since QCD is obtained from QED by extending U(l) gauge group to SU(3) gauge group 
we expect that the factorization theorem in QCD should be valid for light-like Wilson line 
but the factorization theorem in QCD should not be valid for non-light-like Wilson line. 


B. Violation of Factorization Theorem in QED With Non-Light-Like Wilson Line 
in General Axial Gauge 


The correlation function of the type < O|'^(a; 2 )'0(xi) |0 > in the background field method 
of QED in general axial gauge is given by eq. (I45j) which is valid for arbitrary background 
held A >l (x). Hence from eq. (j45l) we find 


< O|'0 (.t 2 )t/’(o:i)|O >a= J [dQ] [dip'] [dip'] {x 2 )i>'{x 1 )e l f d4x ^ 


^DM < ^ M ) 2 +b , [ i 7' ^ 9M-" l + e 7 ,i (^+Q)M]b , ] 


(83) 


because the change of integration variable from unprimed to primed variable does not change 
the value of the integration. 

By using eqs. (17UD and (1771) in (153|) we End 

< 0\$(x 2 )e- ie fo° dXv A ( X2+vX ^e ie fo° dXv ' A ( Xl+vX ^ iIj(xi)\ 0 > a = j[dQ][d#^\ ^(x 2 )^( Xl ) 

Xe i / d4 d-3C?dQ]-^DMQD 2 +#7^- m +e7' i ( J 4+Q), J -e7' i d,4^ v-A(x)]]^ /g^\ 
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Hence from eqs. (178]) . (184]) and (1431) we find 


< 0|^(a; 2 )^(xi)|0 >±< 0\'ijj(x2)e~ ie -fo° dXv - A ( x 2+ vX ) e ie dXv ~ A ( xl+vX )> A (85) 

which implies the violation of the factorization theorem for soft-collincar divergences in QED 
at all order in coupling constant when the non-light-like Wilson line is used in general axial 
gauge. 

C. Violation of Factorization Theorem in QED With Non-Light-Like Wilson Line 
in Light Cone Gauge 

Since eq. (185]) is valid for general axial gauges [see eq. (142]) ] it is also valid for light cone 
gauge (tf = 0). Hence from eq. (185|) we find that in the light cone gauge 

< 0|^(a: 2 )r/;(x 1 )|0 >+< 0\$(x 2 )e~ ie f°°° d ^ A (^+vX) e ie dXv-A^+vX)^^ >a ( g6) 

which implies the violation of the factorization theorem for soft-collincar divergences in QED 
at all order in coupling constant when the non-light-like Wilson line is used in light cone 
gauge. 


D. Violation of Factorization Theorem in QED With Non-Light-Like Wilson Line 
in General Non-Covariant Gauge 


The correlation function of the type < O|'0(x 2 )^(xi)|O > in the background field method 
of QED in general non-covariant gauge is given by eq. (158|) which is valid for arbitrary 
background held A^(x). Hence from eq. (1581) we find 

< 0|^(x 2 )^(a;i)|0 > A = J[dQ][dij/][d^/] i’ , (x 2 )^'(x 1 ) 

e * J d 4 x[-\F% u [Q]-±(2^-d II ,Q„) 2 +4>'[i'y‘ J 'd ll ,-m+eyf 1 (A+Q) tI ,]il>'] /g^\ 


because the change of integration variable from unprimed to primed variable does not change 
the value of the integration. 

By using eqs. (1761) and (1771) in (l87j) we find 

< O|'0(x 2 )e _ d^v-A(x2+ v tygi e / 0 °° |q > ^ = f [^] [^] ${x 2 )^ Xl ) 

xe * I dix \~\ F lv\Q]~4si( I1 ^ d vQA 2 +^[ii^diM-™+ e ‘y ,1 '( A +Q)»- e 7 ,J 'dA^ta V -M X )W\ /gg\ 
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Hence from eqs. (T78|) . (j88|) and (1561) we find 

pOO p oo 

< O|'0(m 2 )'0(^l)|O >^< 0|^(x 2 )e- ie Jo d\vA( X2+ vA) e ieJ 0 dAv.A( X1 +vA)^ Xi ^ 0 >A (gg) 

which implies the violation of the factorization theorem for soft-collinear divergences in 
QED at all order in coupling constant when the non-light-like Wilson line is used in general 
non-covariant gauge. 

E. Violation of Factorization Theorem in QED With Non-Light-Like Wilson Line 
in General Coulomb Gauge 

The correlation function of the type < O|'0(x2)tH x i) |0 > in the background field method 
of QED in general Coulomb gauge is given by eq. (1681) which is valid for arbitrary background 
field A ,l (x). 

Hence by replacing —> [g^ ] everywhere in the derivations in the previous 

sub-section we find that in general Coulomb gauge 

< > 7 ^< 01"0(x 2 ) e~ ie fo° dXv ' A ( X2+vX ) e ie /o°° d\v-A(x 1 +vX ),0 ( Xl ) | q (90) 

which implies the violation of the factorization theorem for soft-collinear divergences in 
QED at all order in coupling constant when the non-light-like Wilson line is used in general 
Coulomb gauge. 

VI. SOFT-COLLINEAR DIVERGENCES AND LIGHT-LIKE WILSON LINE IN 
QCD 

The gauge transformation of the quark field in QCD is given by 

= e i * raua Wip(x). (91) 

Hence one finds that the issue of gauge invariance and factorization of soft-collinear diver¬ 
gences in QCD can be simultaneously explained if cv a (x) can be related to the soft-collinear 
gluon field A^ a (x). This is easily done by using Eikonal Feynman rules in QCD for soft- 
collinear divergences when light-like Wilson line is used which can be seen as follows. 
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The Eikonal propagator times the Eikonal vertex for a gluon with four momentum AT 
interacting with a quark moving with four momentum jT in the limit A <<< p is given by 


15| 


gT a 




(92) 


p ■ A + ie 

In QCD the soft divergence arises only from the emission of a gluon for which all components 
of the four-momentum A M are small (A —>■ 0) which is evident from eq. (l92lh From eq. (1921) 
one also finds that when 0 < A <<< p and p is parallel to A one may find collinear divergence 
in QCD. 

However, since 


_> / vn A 

p-k = p 0 ko - p- k = \p\\k\(J 1 + - cos 6) 


(93) 


one Ends that the collinear divergence does not appear when quark interacts with a collinear 
gluon because of the non-vanishing mass of the quark, i. e., m 7 ^ 0 even if the mass of the 
light quark is very small. From eq. (193|) one Ends that the collinear divergence appears only 
when iji = 0 and 9 = 0 where 6 is the angle between p and A. Since gluons are massless 
and the massless gluons interact with each other one finds that the collinear divergence 
appears in QCD. Since a massless particle is always light-like one finds that the soft-collinear 
divergences can be described by light-like Wilson line in QCD. 

For light-like quark we find from eq. 


gT a 


pH 


= gT a 


l ' 1 


(94) 


p ■ A + ie l ■ A + ie 

where is the light-like four-velocity (|/|=1) of the quark. Note that when we say the ”light¬ 
like quark” we mean the quark that is traveling at its highest speed which is arbitrarily close 
to the speed of light (|/| ~ 1 ) as it can not travel exactly at speed of light because it has 
finite mass even if the mass of the light quark is very small. On the other hand a massless 
gluon is light-like and it always remains light-like. Hence we find that if is light-like four 
velocity then the soft-collinear divergence in QCD can be described by the Eikonal Feynman 
rule as given by eq. (|94l) . Note that the Eikonal Feynman rule in eq. (j94l) is also valid if we 
replace the light-like quark by light-like gluon provided we replace = —if abc . 

From eq. (t94|) we find 


gT a 


d A k l ■ A a (k) 
(27r) 4 l ■ A + ie 


= —gT a i 


dX 


d A k 


e*“ A i • A a (k) = igT a / d\l ■ A"IIX) 


( 95 ) 
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where the gluon field A m (x) and its Fourier transform A^ a (k) are related by 


A» a (x) = 


d 4 k 

Wr 


A^(k)e 


ik-x 


( 96 ) 


Note that a path ordering in QCD is required which can be seen as follows, see also 25]. 
The Eikonal contribution for the soft-collinear divergence in QCD arising from a single 
soft-collinear gluon exchange in Feynman diagram is given by eq. (19511 . Now consider the 
corresponding Feynman diagram for the soft-collinear divergences in QCD due to exchange 
of two soft-collinear gluons of four-momenta /d/ and k 2 . The corresponding Eikonal contri¬ 
bution due to two soft-collinear gluons exchange is analogously given by 

2 r d 4 k x d 4 k 2 T a l ■ A a (k 2 )T b l ■ A b (k l ) 


9 


J (27r) 4 (27 t) 4 (l ■ (ki + k 2 ) + ie)(l ■ h + ie) 

POO POO 

g 2 i 2 / d\ 2 / d\{T a l ■ A a (l\ 2 )T b l ■ A b [l AQ 

Jo J a 2 

r oo poo 

y -V / d\ 2 / dX\T a l ■ A a (lX 2 )T b l ■ A b (l AQ 

Jo Jo 


(97) 


2! jo Jo 

where V is the path ordering. Extending this calculation up to infinite number of soft- 
collinear gluons we find that the Eikonal contribution for the soft-collinear divergences due 
to soft-collinear gluons exchange with the light-like quark in QCD is given by the path 
ordered exponential 

POO 

V exp [ig dXl ■ A a (lX)T a } (98) 

where is the light-like four velocity of the quark. The Wilson line in QCD is given by 


Ve 


ig f* f dx^A^(x)T a 


which is the solution of the equation 26] 


with initial condition 


d^S(x) = igT a A a (x)S(x ) 


S(xi) = 1 . 


When A^ a (x) = A^ a (lX) as in eq. (l98jl we find from eq. 
in QCD for soft-collinear divergences is given by |l3] 


(99) 


( 100 ) 


( 101 ) 


that the light-like Wilson line 


Ve 


ig J x f d2 # M“(a;)T a 


jy e -ig f 0 °° d\l-A a (x f +l\)T a ~\ jy e ig d\l-A b { Xi +lX)T b 


( 102 ) 
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A light-like quark traveling with light-like four-velocity V produces SU(3) pure gauge 
potential A fia (x) at all the time-space position x 9 except at the position x perpendicular to 
the direction of motion of the quark (l ■ x = 0) at the time of closest approach [2-4]. When 
A^ a (x) = as in eq. (195|) we find l ■ x — XI -l — A ^ 0 which implies that the light-like 


Wilson line finds the gluon field A 9a (x) in eq. 
pure gauge is given by 


as the SU(3) pure gauge. The SU(3) 


T a A“(x) = -[d,U{x)} U-\x), 


which gives 


ig 


U(x f ) = Ve* 9 -C dx ' lA ^ x)Ta 


U(x) = e iaTa “ a{ - x) 


U( Xi ) = e i9Taua ( x f\ 


(103) 


(104) 


Hence when A fia (x) = A 9a (Xl) as in eq. (1981) we find from eqs. (11021) and (11041) that the 
light-like Wilson line in QCD for soft-collinear divergences is given by 


Ve 


ie £f dx»A°(x)T a = igT a u{x s ) -igT b A(xi) = 


p e -i9 f 0 °° d\l-A a (x f +l\)T°-~\ p e ig d\vA b { Xi +l\)T 


(105) 


which depends only on end points x 9 and x 9 but is independent of the path. The path 
independence can also be found from the non-abelian Stokes theorem which can be seen as 
follows. The SU(3) pure gauge in eq. (I103p gives 


F‘\A\ = 9„Ai(x) - d„A‘Ax) + gf°*A‘(x)Ai(x) = 0. 


(106) 


Using eq. (11061) in the non-abelian Stokes theorem [27j we find 

Ve i 9 § c dx»A«(x)T“ = Vexp [igj^dx^dx" -p e i9 fy dx ' XAbx(x ' )Tb j F;i(x)T a Ve i9 l 


v dx" 5 A c s (x")T c 


= 1 


(107) 


where C is a closed path and S is the surface enclosing C. Now considering two different 
paths L and M with common end points x f and x^j- we find from eq. (11071) 

Ve ig§ c dx »A°(x)T a = p exp [ig J^Al(x)T a - ig J M dx»A“(x)T a } 


\Ve ig fL 


dxVA“{x)T a 


v e - i9 L dxVAb AVT b 


= 1 


(108) 


which implies that the light-like Wilson line in QCD 


Ve 


ig J Xf dx»A“(x)T a 


(109) 
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depends only on the end points x 9 and x'j but is independent of the path which can also be 
seen from eq. (llOSp . Hence from eq. (1 105 we find that the non-abelian phase or the gauge 
link in QCD is given by 


<£(3.) _ -p e -*9/ 0 °° d\l-A a (x+l\)T a _ e igT a u a (x)' 


( 110 ) 


Note that from eq. (11061) we find the vanishing physical gauge invariant field strength 
square F^ ua [A]F“ w [A] when /D a (x) is the SU(3) pure gauge as given by eq. (11031) . Hence in 
classical mechanics the SU(3) pure gauge potential does not have an effect on color charged 
particle and one expects the effect of exchange of soft-collinear gluons to simply vanish. 
However, in quantum mechanics the situation is a little more complicated, because the 
gauge potential does have an effect on color charged particle even if it is SU(3) pure gauge 


potential and 
simply vanish 


lence one should not expect the effect of exchange of soft-collinear gluons to 
This can be verified by studying the non-perturbative correlation function 
of the type < 0|^(x)^(x / )^(x")'i/’(x w )...|0 > in QCD in the presence of SU(3) pure gauge 
background field. 

Under non-abelian gauge transformation given by 


T a A’“(x) = U(x)T a A a (x)U-\x ) + -[d,U{x)]U~\x), 


W 


U(x) = e i9Tauja{x) (111) 


the Wilson line in QCD transforms as 28] 


Ve 


ie f Xf dx^A' x (x)T a TT/ N 

Jx i M = U[Xf) 


Ve 


ie f* f dx 11 A^(x)T a 


U-\x,). 


( 112 ) 


From eqs. (11051) and (11121) we find 

Ve ~ i9 /o°° dAl-A'“(*+JA)T« = JJ(x)Ve ~ i9 / 0 °° dM-A^x+lX)^, 

which gives from eq. (1 11 01) 


U(x) = exp \igT a u a (x)} 


(113) 


4'(,t) = f/(i)t(,x), 


^(x) = $\x)U 1 (x). 


(114) 


To summarize this, we find that the soft-collinear divergences in the perturbative Feyn¬ 
man diagrams due to soft-collinear gluons interaction with the light-like Wilson line in QCD 
is given by the path ordered exponential in eq. (198|) which is nothing but the non-abelian 
phase or the gauge link in QCD as given by eq. (11 10|) where the gluon field A 9a (x) is 
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the SU(3) pure gauge, see eqs. (11031) . (11041) . (jl05j) . This implies that the effect of soft- 
colliuear gluons interaction between the partons and the light-like Wilson line in QCD can 
be studied by putting the partons in the SU(3) pure gauge background field. Hence we 
find that the soft-collinear behavior of the non-perturbative correlation function of the type 
< 0\^(x)'ijj(x , )'ip(x")'i/j(x" , )...\0 > in QCD due to the presence of light-like Wilson line in 
QCD can be studied by using the path integral method of the QCD in the presence of SU(3) 
pure gauge background held. 

It can be mentioned here that in soft collinear effective theory (SCET) 29] it is also 
necessary to use the idea of background fields B to g.ve well defined meaning to several 
distinct gluon holds [101 ]. 

Note that a massive color source traveling at speed much less than speed of light can not 
produce SU(3) pure gauge held 2-4], Hence when one replaces light-like Wilson line with 
massive Wilson line one expects the factorization of soft/infrared divergences to break down. 


This is in conformation with the hireling in 


% .8] which used the diagrammatic method of 


QCD to prove factorization theorem in NRQCD heavy quarkonium production in case of 


light-like Wilson line [7] and violation of factorization theorem in NRQCD heavy quarkonium 
production in case of non-light-like Wilson line js]. In case of massive Wilson line in QCD 


the color transfer occurs and the factorization breaks down 


8]. Note that in case of massive 


Wilson line there is no collinear divergences which is explained in eq. 


VII. PROOF OF FACTORIZATION OF SOFT-COLLINEAR DIVERGENCES IN 
QCD WITH LIGHT-LIKE WILSON LINE 

By using path integral formulation of QED we have proved in section IV that the factor¬ 
ization theorem is valid in QED when light-like Wilson line is used. Similarly, by using path 
integral formulation of QED we have proved in section V that the factorization theorem 
is not valid in QED when non-light-like Wilson line is used. Since QCD is obtained from 
QED by extending U(l) gauge group to SU(3) gauge group we expect that the factorization 
theorem in QCD should be valid for light-like Wilson line but the factorization theorem in 
QCD should not be valid for non-light-like Wilson line. 

In this section we will prove that the Factorization theorem is valid in QCD for light-like 
Wilson line and in the next section we will prove that the factorization theorem is not valid 
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in QCD for non-light-like Wilson line. 


A. Proof of Factorization Theorem in QCD With Light-Like Wilson Line in Co¬ 
variant Gauge 

The generating functional in the path integral method of QCD in the covariant gauge is 


given by 


30 


Bll 


Z[J,V,v\ = 


A rf ( ^ d ^ la \ J f d 4 x[-±F a ^[Q]-±(d„Qn 2 +W'Y IJ 'd„-m+gT a 7 »QZJ 1 p+J-Q+^+ip v \ 

[ 5(J b J 


where J Ma (x) is the external source for the quantum gluon held Q^ a (x) and r}i(x) is the 
external source for the Dirac held ^(x) of the quark and 


KAQ] = d„Ql(x) - dMx) + gf° tc CfJx)Qi(x), 


flU 


[Q} = F^ a [Q\F a \Q\-( HO) 


determinant det ( ^ ) can be expressed in terms of 


31,1] but we will directly work with the Faddeev-Popov 


Note that the Faddeev-Popov (F-P 
path integral over the ghost helds 
(F-P) determinant det( ^y|jr - ) in this paper. The non-perturbative correlation function of 
the type < O|'0(x2)'0(^i)|O > in QCD in the covariant gauge is given by 


< 0\'ijj(x 2 )^{x 1 )\0 >= J[dQ][d^\[d^\ tp(x 2 )'ip(x 1 ) 


xdet( 


5u b 


) e i f rf4a 


,{Q\-- h (d^Qn 2 -m+gT a ^ Q“] y] 


(117) 


We have seen in the previous section that the soft-collinear behavior of the non- 
perturbative correlation function due to the presence of light-like Wilson line in QCD can 
be studied by using the path integral method of the QCD in the presence of SU(3) pure 
gauge background held. Hence we use the path integral formulation of the background held 
method of QCD in the presence of SU(3) pure gauge background held as given by eq. (I103p 
to prove factorization theorem in QCD. 


Background held method of QCD was originally formulated by ’t Hooft 32] and later 


extended by Klueberg-Stern and Zuber 


33, 


34] and by Abbott 30]. This is an elegant for¬ 


malism which can be useful to construct gauge invariant non-perturbative green’s functions 


in QCD. This formalism is also useful to study quark and gluon production from c 


chromo held 


35] via Schwinger mechanism 


, to compute f3 function in QCD 


assical 


37), to 


(115) 


28 























perform calculations in lattice gauge theories [38] and to study evolution of QCD coupling 


constant in the presence of chromofield 


39]. 


The generating functional in the path integral formulation of the background held method 


of QCD in the covariant gauge is given by 


30 


32j, 


33] 


Z[A, J, r), rj\ = J [dQ} [di/j] [dip] det( 


SG a (Q) 

5u> b 


J J d i x[-^F a l„[A+Q}-^(G a (Q)) 2 +^d l ,-m+gT a ^(A+Qr^+J-Q+f)^+H 


(118) 


where the gauge fixing term is given by 


G°(<3) = d„CT + gK'MF = D^A}CT 


(119) 


which depends on the background held A^ a (x) and 


KM + Q] = 8,[Al + «] - SMI + Ql\ + 9 /*K + Q b J[A: + oil. 


( 120 ) 


We have followed the notations of 


30 


32 


33] and accordingly we have denoted the quantum 


gluon held by Q^ a and the background held by A^ a . Note that in the absence of the external 
sources the SU(3) pure gauge can be gauged away from the generating functional in the 
background held method of QCD. However, in the presence of the external sources the 
SU(3) pure gauge can not be gauged away from the generating functional in the background 
held method of QCD. 

The non-perturbative correlation function of the type < Q\ip{x 2 )ip{xi)\Q >a hi the back¬ 
ground held method of QCD in the covariant gauge is given by 


< O|'0(x 2 )'0(^i)|O >a= J[dQ][dip][dip] ip(x 2 )ip(x 1 ) 


x 


dct( ^ a ^^ )C f d4x l-k Fa ^i A +QG^<,G a m 2 +^d^-m+gT-yS(A + Q)^}^ 


The gauge fixing term A.(G a (Q)) 2 in eq. (111811 [where G a (Q ) is given by eq. (11191) 
invariant for gauge transformation of 


is 


SA a u = gf abc A b u co c + 


(type 1 transformation) 


( 122 ) 


provided one also performs a homogeneous transformation of 30, [33]: 


SQL = gf^QM. 


(123) 
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The gauge transformation of background field A a as given by eq. (11221) along with the 
homogeneous transformation of in eq. (11231) gives 

5(A; + Q“) = gf abc (Al + Q> c + d,co a (124) 

which leaves — \F a2 ^ u [A + Q\ invariant in eq. (11181) . 

For fixed A a , i.e., for 


5A “ = 0 , (type II transformation) 


(125) 


the gauge transformation of 


30 


33]: 


5Q a =gr bc (A b +Q b )u c + d,u a 


(126) 


gives eq. (11241) which leaves — \F a2 [A + Q] invariant in eq. (11181) . Hence whether we use 
type I transformation [eqs. (11221) and (j!23j) ] or type II transformation [eqs. (1125(1 and (11261) ] 
we will obtain the same equation (1 1561) . 

It is useful to remember that, unlike QED js| , finding an exact relation between the gen¬ 
erating functional Z[J, 77, 77 ] in QCD in eq. (11151) and the generating functional Z[A, J,rj,fj\ 
in the background field method of QCD in eq. (11181) in the presence of SU(3) pure gauge 
background field is not easy. The main difficulty is due to the gauge fixing terms which are 
different in both the cases. While the Lorentz (covariant) gauge fixing term — 7 ^(d fl Q ,J ’ a ) 2 in 
eq. (11151) in QCD is independent of the background field A tm (x), the background field gauge 


fixing term — ^-(G a (Q)) 2 in eq. (I118|) in the background field met 


the background field A^ a (x) where G a (Q) is given by eq. (II 19(1 30 


rod of QCD depends on 


32 


. Hence in order 


to study non-perturbative correlation function in the background field method of QCD in 
the presence of SU(3) pure gauge background field we proceed as follows. 

By changing the integration variable Q —> Q — A in the right hand side of eq. (1 121 D we 
find 


< 0\^(x 2 )ip(x 1 )\0 >A= I[dQ][di)\[di)\ 

,dct( ^ ) e* / 

5uj b 

where from eq. (11191) we find 

GJW) - t)/r‘ + 'i.r , 'A'; l cr : - - a,.v“ 


(127) 


(128) 
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and from eq. (11231) [by using eq. (11221) ] we find 




( 129 ) 


The eqs. (11271) . (11281) and (I129|> can also be derived by using type II transformation which 
can be seen as follows. By changing Q —> Q — A in eq. (1121ft we find eq. (11271) where the 
gauge fixing term from eq. (II 19|) becomes eq. (11281) and eq. (11261) [by using eq. (11251) ] 
becomes eq. (11291) . Hence we obtain eqs. (11271) . (11281) and (11291) whether we use the type 
I transformation or type II transformation. Hence we find that we will obtain the same eq. 
(11561) whether we use the type I transformation or type II transformation. 

Changing the integration variable from unprimed variable to primed variable we find from 


eq. (H 27D 



This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

The equation 


q';(x) = + gf^(x)Qi( X ) + a^(x) 


(131) 


in eq. (II 29ft is valid for infinitesimal transformation (oj « 1) which is obtained from the 
finite equation 


T°Q'Z(x) = U(x)T“QZ(x)U-‘(x) + -\d,U(x)}U-\x) 

*9 


(132) 


where 


U(x) = 'p e -i 9 f 0 °°d\l-A^x+l\)T“ 


= e igT°-uj°-{x) 


(133) 


Simplifying infinite numbers of non-commuting terms we find 



(134) 


where 


M ab (x ) = f abc u c (x). 


(135) 
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Hence by simplifying infinite numbers of non-commuting terms in eq. (11321) [by using eq. 
(1134p and j^]] we find that 

p gM(x) _ 1 

Q'» = + [ • U [S^\x)\, M ab {x) = r^(x). (136) 

Under the finite transformation, using eq. (I136j) . we hnd 

IdQ'} = [dQ} det[^] = [dQ\ det[[e° M ^}} = [dg]exp[Tr(ln[e^ M W])] = [dQ\ (137) 
where we have used (for any matrix H ) 

dettf = exp[Tr(lntf)]. (138) 


Let us define 


y/(x) = Pe- i9 /»” JA '- A “<“+ IA > T >(x) = e‘ 9T "““ (i y(x). 


(139) 


When A^^x) is SU(3) pure gauge we hnd by using eqs. (11361) and (I139|) that 


[#'][<#'] = [dii][dip], V>'[n'‘d u -m + sTYQZW = - m + 


KAQ'} = F a lM 


(140) 


Simplifying all the inhnite number of non-commuting terms in eq. (I103[) we hnd that the 
SU(3) pure gauge A ,M (x) is given by j^] 

' e gM(x) _ Y 


A^ a (x) = d»u b ( x) 


gM(x) 


ab 


where M ab (x ) is given by eq. (11351) . From eqs. (11411) and (I136p we hnd 


(141) 


Q'» - K(x) = [ e 9M(x) UQ*(z), M^x) = r*V’(i). 


(142) 


Using eqs. (I137p . (I140p and (I142p in eq. (I130D we hnd 

< 0|^(x 2 )^(^i)|0 >a= I[dQ}[d^][di^(x 2 )[Ve~ l9T h ^^^(x i) 


x 


det (f3^1) e i J d 4 xl-jF^[Q]-^(G-(Q')) 2 +^[i 7 ^-m+gT^Q-] 


6u b 

From eq. (I128P we hnd 


(143) 


GW) = 4//'“ + gf^A'-d - c(.4'“. 


(144) 
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By using eqs. (11361) and (1141)) in eq. (1144)) we find 


p gM(x) _ -| 

G a f{Q') = dV M ^] ab Q b Jx) + [ U [d,u\x)]] 


+gf abc [d fl uj e (x) 
-d^u) b (x 
which gives 


e gM(x ) _ 

gM ( x ) 

e gM{x) _ ^ 


gM(x ) 

P gM(x) _ i 

+ 1 g M ( X ) 1 ^ 1 ^ MU 


gM ( x) 


( 145 ) 


- ab 


G° f (Q') = SV^UQ^x)] 

e gM(x) _ ^ 


+gf° bc {dV(x) 

From eq. (11461) we find 


gM (x) 


gM(x) _ -| 

\Ue” MM UQt(x) + [ ' gM(x) U [V*W]]- (146) 


G%Q') = d»{{e’ M <%„Ql(x)} + gr^u’ix)] 

which gives 


e gM(x) _ 2 

gM(x) 


]{{e‘ M ^UQ d Jx)](U 7 ) 


be 


G a f (Q') = [e sM( % b d^(x) 

+ Q b ^ x )dV M ^} ab } + \[d^\x)} 

From we find 


e gM(x) _ ^ 

gM (x) 


] gf*'[[ e *WI] al QZ(x)l. 


(148) 


gV^igTVfz)].. = f^fl)] 


e gM(x ) _ ^ 

gM ( x ) 


35 ^%, M ab (x) = r bc u c (x) 


ab 


(149) 


which in the adjoint representation of SU(3) gives (by using T fe “ = —if abc ) 

' e gM{x) _ 1 ' 


^e 9M ( x )] ad = [d^oj e (x)\ 


gM ( x) 


gf bac [e M ^] cd , Mab(x) = f abc co c (x). (150) 


be 


Using eq. (I150p in (I148p we find 


GVQ 1 ) = [ ( f M ^Ud"Q t Jx) 


(151) 


which gives 


(' G a f m 2 = {d,Q^{x)f. 


(152) 
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Since for n x n matrices A and B we have 


det(AB) = (detH)(det£>) 


(153) 


we find by using eq. (j 151 j) that 

det[^)] = = det V^]J (d ^ X)) \ 


= det[ 


det[[e^U 

5(d^ a (x)) 


det [ 


5(d^Q c ^(x)) 

5u b J 


= exp[Tr(ln[e ffM ^ x) ])] det[ 


5(d,Q^(x)). 

5uj b 


Suj b 


(154) 


Using eqs. (II52ft and (11541) in eq. (1143ft we find 

< O|'^(x 2 )'0(zi)|O >A= J[dQ][d^\[d^\ i>(x 2 )[Pe~ 1 ' 9T h 


x 


d Ct f ^(^^°) ) e i f ^A-i^JQl-^idvQ^+AB^-rn+gT^QZM 
V S0J b ' 


(155) 


Using the similar technique as above we find 

< 6\^)(x 2 )[Pe 9T h dTMj4,l(x) ]^(xi)|0 > A = j[dQ][dip][dip] ^(x 2 )^(x 1 ) 


x 


dct( ^(^ a ) ) e i f d4 A- 2 F a lAQ}-^(dvQn 2 +j[B»d^-m+gT"-^Q“} 
v 5u b 


(156) 


in the presence of SU(3) pure gauge background field A^ a (x) as given by eq. (1103(1 where 
M ab (x) is given by eq. (11351) . 

Note that eq. (1 156ft is valid whether we use type I transformation [eqs. (1122(1 and (112311 
or type II transformation [eqs. (1125ft and (1126(1 ] . However, since eq. (II 1 1(1 is used to study 
the gauge transformation of the Wilson line in QCD as given by eq. (11121) . we will use 
type I transformation [see eqs. (11221) and (11231) ] in the rest of the paper which for the finite 
transformation give eq. (11111) and [301. [33] 

T a Q'“(x) = U(x)T a Ql(x)U- 1 (x), U(x) = e* 9W(a;) . ( 15 7) 

From eq. (11571) we find that under gauge transformation given by eq. (II 1 1(1 the (quantum) 
gluon field Q fia (x) transforms as 


«?(*) = l^UQ^x) 

where M ab (x) is given by eq. (11351) . From eqs. (11171) and (1156(1 we find 

< 0|^(x 2 )V'(^i)|0 >=< Q\'ip{x 2 )'(Pe gT <i ^' 4 ' ,(:rl ]V'(xi)|0 > A 


(158) 


(159) 
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which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD when light-like Wilson line is used in the covariant gauge where 


igT a f^dx»A“{x) = [ Vexp [_ ig r dXl . A C( X2 + 

JO 


Ve 

roo 


xVexp[ig J dXl ■ A b (x± + lX)T b ] 


(160) 


is the light-like Wilson line in terms of the non-abelian gauge link or non-abelian phase and 
is the light-like four velocity. 


B. Proof of Factorization Theorem in QCD With Light-Like Wilson Line in Gen¬ 
eral Axial Gauge 


In QCD the generating functional with general axial gauge fixing is given by 20, [21] 


Z[J, V, rj\ = 


xe 


i J d i x[-^F a l„[Q]-^ ; (ri l _ l Q> la ) 2 +ip[i'y' J -d l _ l -m+gT a 'yl J 'Q^]ip+J-Q+'il)ri+rjip] 


(161) 


where is an arbitrary but constant four vector 


? 7 2 < 0, axial gauge 

rj 2 = 0, light — cone gauge 

7] 2 > 0, temporal gauge. 


(162) 


Note that unlike covariant gauge in eq. (II 15jl there is no Faddeev-Popov (F-P) determinant 


211. The non¬ 


in eq. (I16ip because the ghost particles decouple in general axial gauges 
perturbative correlation function of the type < O|'0(x2)'*/ ; (^i)|O > in QCD in general axial 
gauge is given by 


< O|'0(rr 2 )'0(^i)|0 >= J[dQ ][#][#] ^{x 2 


xe 


i f d 4 x[-±F a2 l/ [Q\-±(r l VQa) 2 +i/,[i 7 vd /J ,-m+gT a 'yf J -Qa 


(163) 


The generating functional in the background field method of QCD with general axial 
gauge fixing is given by [21] 


Z[A, J,rj,rj\ = 


xe 


i f d 4 a;[-iF a ^[A-|-Q]-^(r)' i Q“) 2 +i/;[i7^9 M -m+ 9 r a 7' J (n+Q)“]V>+J-Q+r/y+i/’d 


(164) 
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The non-perturbative correlation function of the type < Q\ip(x 2 )ip(x\)\Q >a in the back¬ 
ground field method of QCD in general axial gauges is given by 


< 0\ip(x 2 )ip(x 1 )\0 > A = J[dQ][dip\[dip] ippx^ippx^) 


xe 




(165) 


By changing the integration variable Q —> Q — A in the right hand side of eq. (I165|) we 
find 


< 0\'if}(x 2 )'tjj{xi)\0 >A= J [dQ][dip][dip] ip(x 2 )ip{x l ) 


xe 


i f rf 4 rr[-| F *2j Q] -£ ( ^ ( Q-A)Z)2+rti 7 v d/x - m+gT °- 7l i Q ° 


(166) 


Changing the integration variable from unprimed variable to primed variable we find from 
eq. (11661) 


< 


O|'0(^ 2 )'0(^i)|0 >n= f [dQ'][dip'][dip'] ip'(x 2 )ip' (aq) 


xe 


i f d 4 x[-±F° 2 ^[Q']-±(r,v(Q'-A)Z) 2 +^li 7 v dM -m+gT« 7 »Q'“W'] 


(167) 


This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

Using eqs. (I137[) , (I140p and (j!42|) in eq. (I167P we find 

< 0|^(x 2 )^(^i)|0 >a= J [dQ][dip][dip] ip(x 2 )[Ve~ l9T h i) 

Xe*/ d4x [-l Fa lM-^W([e 9M ( x PabQ b ^(x))) 2 +$[i7»d»-rn+gT a 7 i 1 -Q“]'il>] (168) 


which gives 


< 0\ip(x 2 )ip(xi)\0 > A = J [dQ] [dip] [dip] ip(x 2 )[Pe igT h dx ^ Al " (x) ]ip(x 2 ) 


i J d 4 x[-±F“l„[Q]-±WQ“) 2 +j\i^d^-m.+gT"-^Q“} 


xe 

Using the similar technique as above we find 

< 0| ip(x 2 )[Ve 9T h dx>1A ^ x) ]ip^ x i) |0 >a= J[dQ][d'ip][dip) '4){x 2 )ip{x 1 ) 

xe i f d 4 x[-\F a l v \Q]-^(^Ql) 2 +^a^- m +gT^Ql}^} 


(169) 


(170) 


in general axial gauge in the presence of SU(3) pure gauge background field /U“(x) as given 
by eq. (|103p . 
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From eqs. (I105j> . (I163j> and (I170H we find 


_ _ /*oo 

< 0\'ijj(x2)'ip(xi)\0 >=< 0\'ip(x2)[Vexp[—ig / dXl ■ A c (x 2 + l\)T c ]\ 

Jo 

poo 

xVexp[ig J dXl ■ ^{xi + /A)T 6 ]'0(£i)|O >a (171) 

which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD when light-like Wilson line is used in general axial gauge. 


C. Proof of Factorization Theorem in QCD With Light-Like Wilson Line in Light 
Cone Gauge 


The light-cone gauge corresponds to [20-22] 


v ■ Q a = o, 


rf = 0 


(172) 


which is already covered by eqs. (116111 and (11621) where the corresponding gauge fixing term 
is given by . In the light-cone coordinate system the light-cone gauge 22] 


Q +a = 0 


(173) 


corresponds to 


V M = i.V + ,V ,V±) = (0,1,0) (174) 

which covers ij ■ Q a = 0 and rj 2 = 0 situation in eq. (11721) . 

Since eq. (1 171 [1 is valid for general axial gauge [see eq. (1162(1 ] it is also valid in light cone 
gauge (rf = 0). Hence in light cone gauge we find from eq. (117111 that 

< 0\'ijj(xi)'i/j(x2)\0 >=< 0\'i/j(xi)[Vexp\—ig / dXl ■ A c (x i + IX)T C ]\ 

Jo 

poo 

xVexp [ig J dXl ■ A b (x 2 + lX)T b ] r i/j(X2)\0 >a (175) 

which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD when light-like Wilson line is used in light-cone gauge. 
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D. Proof of Factorization Theorem in QCD With Light-Like Wilson Line in Gen¬ 


eral Non-Covariant Gauge 


In QCD the generating functional with general non-covariant gauge fixing is given by 


23 


24] 


Z[J,V,v\ = 


S(^d,Q 

det(- 1 M 


8uj b 


xe 


i f d 4 xl-jF a2 u JQ]-^(^fd^) 2 +^li 7 ^-m+gT^Q^+J-Q+^rj+ft^ 


(176) 


where is an arbitrary but constant four vector. The non-perturbative correlation function 
of the type < |0 > in QCD in general non-covariant gauge is given by 


< 0\'G(x 2 )'ip(x 1 )\0 >= J[dQ][di>}[di>\ ip(x 2 )i/;(x i) 


x 


det(», D-H- 


- 3 F a lAQ] -h (■ dn-Ql) 2 -m+gT^YQ“ 


(177) 


The generating functional in the background held method of QCD with general non¬ 
covariant gauge fixing is given by 23,124] 

,5G a (Q), 


Z[A,J,rj, rj\ = J[dQ\ [dip] [dip] det(- 


5u b 


-) 


xe 


i J d 4 x[-iF“2 i/[ A+Q ] -i ( e“ ( Q) ) 2+^[i 7 Ma^_ m+g T“ 7 /'(A+Q ) “ ] V'+J-Q+^+^] 


where 


S“(Q) = Cf + gf^KQi) = -4r D A A \Ql 


(178) 


(179) 


r) " r)- 

The non-perturbative correlation function of the type < 0\ip(x 2 )ip(xi)\0 >a in the back¬ 
ground held method of QCD in general non-covariant gauge is given by 


< 0\ip(x 2 )‘ip(x 1 )\Q > a = 


ip(x 2 )ip(x l) 


X 


dct( ^f )C f d ^i-\ Fa lA A +Q] -^(G a m 2 +^d^-m+g T “ 7 ^A+Q)-W]^ ( 18 Q) 


By changing the integration variable Q —> Q — A in the right hand side of eq. (I180|) we 
hnd 


< 0\ip(x 2 )ip(x 1 )\0 >a= J [dQ][dip][dip] ip(x 2 )ip(x 1 ) 


x 


dct( ^ ) e i / d4;r [-3 F “^[<3]-^(^(Q)) 2 +V ; [n' 1 ^-m+gT“ 7 MQ“]v,] 


(181) 
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where from eq. (I179j> we find 


Gf(Q) = 




(d,Ql + gr c KQv - - -2 9f abc (v ■ A b )(v ■ A c ) 


” mAm-^a.A: 


( 182 ) 


g* 7 

Changing the integration variable from unprimed variable to primed variable we find from 
eq. (1 181 j) 

< 0 \^{x 2 )i>{x 1 )\0 > A = J [dQ'}[di’'][d^'\ ^(x 2 W{x l ) 

xdet(^^-^) e * f ( 183 \ 
y 5u b ' 


This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

Using eqs. (1 13 71) , (11401) and (I142j) in eq. (I183j) we find 


< 0|^(x 2 )V'(^i)|0 >a= J[dQ][di/j\[difi\ ip(x 2 )[Ve igT 


x 


[ Scu b J 


(184) 


From eq. (11821) we find 


rfrf 


etc?) = 2 - 2 -[a„g? + gf^A^Q':] - "M. a„Ai. 


By using eqs. (11361) and (1 141 p in eq. (j!85p we find 


g v p gM(x ) _ i 

etm = ^m^UQKz) + I—777^-U [d„?(x)}} 


r y 


+gr bc [d,Lu e (x) 


0 gM(x) _ ^ 

gM (x) 


be 


gM (x) 

p gM(x) _ 1 

[e lMW ] c d Qt(x) + \ x) U 


rj^rj 1 ' 


d^[d v u b (x) 


0 gM(x) 


gM (x) 


ab 


which gives 


(185) 


(186) 


e,“W') = ^[a,[[e»"<*>UQ‘(*)] 


+gf" bc [d„u’(x) 


e gM(x) _ ^ 


gM(x) 


- be 


p gM(x ) i 

{e« Ml % d Qt(x) + [ g U [0^(x)]]]. (187) 
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From eq. (I187j> we find 


Sf(Q') = ^UQtw] 




e gM(x) _ ^ 

gM{x) 


W’ M{x, UQi(x)\A 


(188) 


- be 


which gives 


«/«') = ^He^UdMx) + Q b „(x)3„{[e‘ M <% b } 

p gM(x) _ i “I 

}gf" bc [[e sm ‘ ) UQi(x)}}. 

be 


+[[d li uj e (x)\ 
Using eq. (I150p in (I189p we hnd 


gM(x) 


(189) 


Sf(Q') = ^ [e 9 " W ]»i AQlM 


which gives 


, rfrj v 


(G a f m = (^:w ) 2 


(190) 


(191) 


From eqs. (11531) and (11901) we hnd 


Hrtr s Gf(Q') ] _ H .\^V" $[[e 9M{x) ]acd li Q c l ,( x )\ ] _ , rJ/Yr.jtff*), Kdg,Ql(x)) 
d6t[ 5oj b ] ~ [ ry 2 5u b 1 ~ [ ? [ Jac Soj b 1 


det[[e : 


gM(x) 


dct ^V u Sjd^QKx)) . 


V 


5uj b 


= exp[Tr(ln[e»"W])] det |!!^W«). 


5uj b 


= dct[ ^^ M x)). 


rj* Sco b 

Using eqs. (1 191 p and (I192p in eq. (I184P we hnd 


(192) 


< 


O|' 0 (a: 2 )^(a:i)|O >, 4 = J[dQ][difi][diJi\ i/j{x 2 )[Pe l9T h ^^^^(^i) 


x 


8 { d nQl) ^ J [ d 4 xl-iF^JQ]-^(^fd / xQ-) 2 +^[ira / x-m+gT-rQZli’ 
deH ^ 5u b )G 

Using the similar technique as above we hnd 

0| 4>{x 2 )[Pe 9T dx ^ A ^ x) ]^(xi)\0 > A = j[dQ ][#][#] $(x 2 )ip(xi) 


(193) 


< 


x 


dct ^V \ e » /d 4 a;[-fF a ^[Q]-^(2hj_9 /1 Q“) 2 +^[j 7 M9 /J -m+9T a 7''Q“]^] n 94) 

1 rf 5oj b { 


in general non-covariant gauge in the presence of SU(3) pure gauge background held A^ a (x) 
as given by eq. (11031) . 
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From eqs. (I105j> . (11771) and (j!94|) we find 

_ _ poo 

< 0\'ijj(x9) r ip(xi)\0 >=< 0 \^(x 2 )[Vexp[—ig / dXl ■ A c (x 2 + l\)T c ]] 

Jo 

roo 

xVexp[ig J d\l ■ A b (xi + l\)T b ]iJj(Xi)\0 >a (195) 

which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD when light-like Wilson line is used in general non-covariant gauges. 


E. Proof of Factorization Theorem in QCD With Light-Like Wilson Line in Gen¬ 
eral Coulomb Gauge 


In QCD the generating functional with general Coulomb gauge fixing is given by 23, 24] 


Z[J, If, If] = 


det(- 


5u b 


xe 




(196) 


where 


nA = ( 1 , 0 , 0 , 0 ). 


(197) 


The non-perturbative correlation function of the type < 0\'if>(x 2 )'tp( < xi)\0 > in QCD in general 
Coulomb gauge is given by 


< O|’0(x 2 )'0(a ; i)|O >= J[dQ][dip][d'tjj] ip(x 2 )ip(x 1 ) 


x det(—- n 


Suj b 


2 }dgQl) ^ e i J d i x[-\F-l v [ Q ]-^([g^-^)d^Ql) 2 +^d^-m+g T a ^ Q -\ 


(198) 


The generating functional in the background field method of QCD with general Coulomb 


gauge fixing is given by 23|, 


Z[A, J, rj, rj\ = 


d t(— 


xe 


i f d 4 x[-±F a ? l JA+Ql-£(g a (Q)) 2 +Tf,[i'y»d /J ,-m+gT a 7 v(A+Q)Z]Tp+J-Q+fj^+TpT ) \ 


(199) 


where 




rrn 


»n v 


G a (Q ) = \gT ~ —rlidM + gf abc A b Ql) = fr 




n 


n^n 


rr 


}d,[a]q: 


( 200 ) 
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The non-perturbative correlation function of the type < >a in the back¬ 

ground field method of QCD in general Coulomb gauge is given by 


< 0|^(x 2 Mah)|0 >j4= J[dQ][di/j] [cb/’] ^{x 2 )^{x i) 


x 


dctf ^^ V f d ^[~\ Fa l^ A +Q]-^ a m 2 +^d»-rn+gT a ^{ A +Q)“]-il>} ^ 2 01) 

S0J b ' 


Hence by replacing ^->■ [g^ — everywhere in the derivations in the previous sub¬ 

section we find 


_ _ r oo 

< 0\'ijj(x 2 ) r ip(xi)\0 >=< 0 \'ip(x 2 )[Vexp[—ig / dXl ■ A c (x 2 + l\)T c ]] 

•/ii 

r oo 

xVexp[ig J dXl ■ A b (xi + |0 (202) 

which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD when light-like Wilson line is used in general Coulomb gauge. 


VIII. VIOLATION OF FACTORIZATION THEOREM OF SOFT-COLLINEAR 
DIVERGENCES IN QCD WITH NON-LIGHT-LIKE WILSON LINE 


As mentioned earlier, in some of the studies, the non-light-like Wilson line is used in 
the definition of the transverse momentum dependent parton distribution function (TMD 
PDF) at high energy colliders jl, 6 ]. These studies involve diagrammatic calculation at one 
loop level using perturbative QCD. However, the (transverse momentum dependent) parton 
distribution function and fragmentation function are non-perturbative quantities in QCD. 
Hence if one uses the perturbation theory to study their properties l], 6[, then one may end 
up finding wrong results. In general, the non-perturbative phenomena may be impossible to 
understand by perturbation theory, regardless of how many orders of perturbation theory 
one uses. Take for example, the non-perturbative function 


f(x) = e 


— c> 


(203) 


The Taylor series at x = 0 for this function f(x ) is exactly zero to all orders in perturbation 
theory, but the function is non-zero if x ^ 0. 

Hence the diagrammatic method using perturbative QCD may not be always sufficient 
to prove factorization of soft-collinear divergences of the (transverse momentum dependent) 
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parton distribution function and fragmentation function at high energy colliders, regardless 
of how many orders of perturbation theory one uses, because the (transverse momentum 
dependent) parton distribution function and fragmentation function are non-perturbative 
quantities in QCD. 

On the other hand the path integral method of QCD can be used to study non- 
perturbative QCD. The path integral method of QCD is valid at all order in coupling 
constant. Hence we find that the path integral method of QCD is very useful to prove 
factorization of soft-collinear divergences of non-perturbative quantities in QCD such as the 
(transverse momentum dependent) parton distribution function and fragmentation function 
at all order in coupling constant at high energy colliders. 

In this paper we have proved in section VII, by using path integral method of QCD, that 
the factorization theorem in QCD is valid at all order in coupling constant when light-like 
Wilson line is used. In this section we will prove, by using path integral method of QCD, 
that the factorization theorem in QCD is not valid at all order in coupling constant when 
non-light-like Wilson line is used. This implies that the factorization theorem is violated 
in all the previous studies jl, 6 | which used the non-light-like Wilson line in the definition 
of the (transverse momentum dependent) parton distribution function and fragmentation 
function at high energy colliders. This is in conformation with the finding in [7], 8] which 
proved factorization theorem in NRQCD heavy quarkonium production in case of light-like 
Wilson line [7| and the violation of factorization theorem in NRQCD heavy quarkonium 
production in case of non-light-like Wilson line js]. In case of massive Wilson line in QCD 
the color transfer occurs and the factorization breaks down jg]. The simple physics reason 
behind this is the following. The light-like Wilson line produces pure gauge potential which 
gives F^Jx) = 0 which can be gauged away in the sense of factorization because pure gauge 
corresponds to unphysical longitudinal polarization and hence factorization theorem works 
(see sections IV and VII for details). The non-light-like Wilson line does not produce pure 
gauge potential (F° M (x) ^ 0), the form of which can be arbitrary which can not be gauged 
away in the sense of factorization because non-pure gauge does not correspond to unphysical 
longitudinal polarization and hence the non-light-like Wilson line spoils the factorization (see 
sections V and this section for details). 
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A. Violation of Factorization Theorem in QCD with Non-Light-Like Wilson Line 


In Covariant Gauge 


Note that eq. (11271) is valid for any arbitrary background field A^ a (x). Hence we start 
from eq. (11271) by following exactly the same procedure that was followed for light-like 
Wilson line case. By changing the integration variable from unprimed variable to primed 
variable we find from eq. (11271) 

< O|-0(rc 2 )'0(^i)|0 >a= J[dQ'][dft][dft] ft(x 2 )ft Oi) 

.dctf *^) ) e*/ d4x l-i Fa UQ'}-£( G f(QA 2 +j'[irdv- m +gT^Q'ZW} (204) 

v 5u b ’ ■ \ ) 

This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

Now let us define the primed variables Q' and ft for the non-light-like Wilson line case 
in analogous to the light-like Wilson line case. First of all in analogous to eq. (I139|) for the 
light-like Wilson line case let us define the primed variable ft for the non-light-like Wilson 
line case 

ft(x) = Ve~ i9 fo°° d\v. A °(x+v\)T«^ (205) 

where A^ a [x) is not the SU(3) pure gauge and v 9 is the non-light-like four velocity. Similarly 
in analogous to eq. (I132|) for the light-like Wilson line case let us define the primed variable 


Q' for the non-light-like Wilson line case 


T a Q’ftx) = U(x)T a Q a ftx)U-\x) + 7 -[d li U{x)]U~ 1 (x) 


1 

m ' 


(206) 


where 


U(x) =p e - i9 fo°° d * v - Aa ( x + v V Ta (207) 

where A fia (x) is not the SU(3) pure gauge and v /J ' is the non-light-like four velocity. Under 
this finite transformation, using eq. (I206P . we find 

yo 

(208) 


[dQ'] = [dQ} det[^7] = [dQ] 


dQ b 


which is similar to eq. (I137jl for the light-like Wilson line case. From eqs. (12061) . (I207p and 
(12051) we find that 

[dft] [dft] = [dft [dft], ft[iftd^ -m + gT a ftQ'ftft = -m + gT a ftQl]ft 


F a LlQ'} = KM 


(209) 
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for the non-light-like Wilson line case which is similar to eq. (j!40|) for the light-like Wilson 
line case. 

Using eqs. (I208P and (12091) in eq. (12041) we find 


exp 


< 0\i>(x 2 )^(xi)\0 >A= 

poo 

xPexp \—ig / dXv ■ A b (x i + v\)T b ]]iJ;(xi) 

Jo 

J I d 4 x{-±F“l„[Q]-±(G a f (QX) 2 +i’[i7»dp-rn + gT^Q“}Tp} 

1 5u b ’ 


dXv ■ A c (x 2 + uA)P c ]] 


X 


From eq. (1128ft we find 

G a f (Q') = d /l Q' ,m + gf abc AlQ^ c - d^ a . 

Now from eqs. (1206ft and (1207ft we find that for the non-light-like Wilson line case 
<3'» - A‘(x) / N ab (x) = /“ *W(x) 


( 210 ) 


( 211 ) 


( 212 ) 


when A^ix) is not the SU(3) pure gauge where h a (x) is any function. Hence we find that 
eq. (1212ft for the non-light-like Wilson line case differs from the corresponding eq. (11421) for 
the light-like Wilson line case. From eqs. (12111) . (12061) . (1207ft and (1212ft we find 


G%Q>) + 

where w a {x ) is any function which gives 

(g?(o '» 2 + (a„<r M) 2 


Sa,(x) = f a w‘(x) 


(213) 


(214) 


for the non-light-like Wilson line case which differs from the corresponding eq. (1152(1 for the 
light-like Wilson line case. From eq. (12131) we find 


det [ 




(215) 


Soj b Su> b 

for the non-light-like Wilson line case which differs from the corresponding eq. (1154(1 for the 
light-like Wilson line case. 

Hence we find from eqs. (12101) . (1214ft and (121511 that 

/»nn 

dXv ■ A c (x 2 + vX)T c }} 


(216) 


< 0\ip(x 2 )^(x 1 )\0 >at^ 

poo 

xPexp \—ig / dXv ■ A b (x i + vX)T b }]i/j(xi) 

Jo 

e i f di A-\F a lAQ]-^ d ^A 2 +^d^-rn+gT^Ql 

v Su b 
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which gives 


< 0\ijj(x2)[[Pexp[—ig / dXv ■ A c {x 2 + r?A)T c l] x 
Jo 


o 


dXv ■ A b (x\ + vA)T 6 ]]'0(xi)|O >_4 


7 ^ j [dQ)[dft[dfti>(x 2 )i>{xi) det( 


Six b 


) e i f d * 


A-\F a lAQ}-h(W l ) 2 +‘4> [ W d„-m+gT°-,»Q“ 


(217) 


for the non-light-like Wilson line case which differs from the corresponding eq. (I156|) for the 
light-like Wilson line case. 

Hence from eqs. (II 17(1 and (1217(1 we hnd 

< 0\^(x 2 )ip(xi)\0 > 

7 ^< 0 \'il>(x 2 )[['Pexp[—ig J dXv ■ H c (x 2 + uA)T c ]] x Pexp [ig J dXv ■ A b (x i + nA)T 6 ]]^(xi )|0 >a 

(218) 


which proves that the factorization theorem is violated at all order in coupling constant when 
the non-light-like Wilson line is used in QCD in covariant gauge where v M is the non-light-like 
four velocity. 


B. Violation of Factorization Theorem in QCD with Non-Light-Like Wilson Line 
In General Axial Gauge 

The non-perturbative correlation function of the type < 0|$(xi)^(a: 2 )|0 in the back¬ 
ground field method of QCD in general axial gauge is given by eq. (I166P which is valid for 
any arbitrary background field A^ft). Changing the integration variable from unprimed 
variable to primed variable we find from eq. (I166D 

< 0\^(x 1 )'ip(x 2 )\0 >A= I [dQ'l [dft] [eh//] ft(xi)ft(x 2 ) 

Xe iJ (219) 

This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

Now from eqs. (I209jl . (j2 12|) . (12 19[) and (I163jl one finds in the general axial gauge 

< 0\ftx 2 )lp(x 1 )\0 > 

_ POO POO 

7 ^< 0 \'il’(x 2 )[['Pexp[— ig J dXv ■ A c (x 2 + vX)T c ]] x Pexp [ig J dXv ■ A b (x i + vX)T b ]]ilj(xi)\0 >a 

( 220 ) 
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which proves that the factorization theorem is violated at all order in coupling constant 
when the non-light-like Wilson line is used in QCD in the general axial gauge where is 
the non-light-like four velocity. 

C. Violation of Factorization Theorem in QCD with Non-Light-Like Wilson Line 
In Light Cone Gauge 

Since eq. (12201) is valid for general axial gauges [see eq. 11162(1 ] it is also valid for light 
cone gauge (r/ 2 = 0). Hence from eq. (12201) we find that in the light cone gauge 

< O|'0(cc 2 )'0(^i)|0 > 

_ poo poo 

7 ^< 0\'i/;(x 2 )[[Vexp[—ig / dXv ■ A c (x 2 + v\)T c ]] x "Pexp [ig / dXv ■ A b (x i + vX)T b ]]tp(xi)\0 

Jo Jo 

( 221 ) 

which proves that the factorization theorem is violated at all order in coupling constant 
when the non-light-like Wilson line is used in QCD in the light cone gauge where v ,J ' is the 
non-light-like four velocity. 


D. Violation of Factorization Theorem in QCD with Non-Light-Like Wilson Line 
In General Non-Covariant Gauge 


The non-perturbative correlation function of the type < O|'0(xi)'0(x 2 ) |0 in the back¬ 
ground field method of QCD in general non-covariant gauges is given by eq. (11811) which 
is valid for any arbitrary background field A^ a (x). Changing the integration variable from 
unprimed variable to primed variable we find from eq. (1181ft 


< O|^(xi)'0(x 2 )|O > A = 


'] i/j , (x l )'ip'(x 2 ) 


X 


dct( ^) ) e i / dM-zF a lAQ']-£;(G]m) 2 +Vlird»-m+gT“rQ'“W] 

1 5u b 1 ’ 


This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. From eqs. (I182p . (1206|) . (1207|) and (1212[) we 
find 


Q‘,(Q') + ( d »QU x )> 


Sat(x) = r c w“(x) (223) 
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where w a (x) is any function which gives 


(S a ,m 2 * WMf (224) 

for the non-light-like Wilson line case which differs from the corresponding eq. (I19ip for the 
light-like Wilson line case. From eq. (I223j) we find 


det[ 


sg a f (Q') 

5u b 


8( 

) * det[^L_ 


d^Q a Ax)) 

5co b J 


(225) 


for the non-light-like Wilson line case which differs from the corresponding eq. (1 192 p for the 
light-like Wilson line case. 

Now from eqs. fl209f) . (I212j) . (12241k (12251k (12221) and (11771) we find that in the general 
non-covariant gauge 


< O|'0(x 2 )'0(^i)|O > 

_ poo poo 

7 ^< 0 \^(x 2 )[[Pexp[—ig / dXv ■ A c (x 2 + v\)T c ]\ x Vexp[ig / d\v ■ A b (x i + v\)T b ]]'ijj( < xi)\0 

jo Jo 

(226) 


which proves that the factorization theorem is violated at all order in coupling constant when 
the non-light-like Wilson line is used in QCD in the general non-covariant gauge where 
is the non-light-like four velocity. 


E. Violation of Factorization Theorem in QCD with Non-Light-Like Wilson Line 
In General Coulomb Gauge 

The non-perturbative correlation function of the type < 0|^(xi)'i/’((C2)|0 >u in the back¬ 
ground field method of QCD in general Coulomb gauge is given by eq. (120 ip which is valid 
for any arbitrary background held A^ a (x) where 

rd l n v . 


nu 


g a {Q) = [<r 


n^n 


rr 


-}{d,Ql + gr c A b Ql) = [g 


PLV 


n* 


}D,[A\Ql- 


(227) 


and 


n l “ = (1,0, 0,0). 


(228) 
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Hence by replacing ^-* [g^ u — ' n ^r] everywhere in the derivations in the previous sub¬ 

section we find that in the general Coulomb gauge 

< 0\^(x 2 )^(x 1 )\0 > 

_ POO POO 

7 ^< 0\'il>(x2)[[Pex'p[—ig / dXv ■ A c {x 2 + v\)T c ]] x Vexp[ig / d\v ■ A b (x i + vX)T b ]]'ip(xi)\0 

Jo Jo 

(229) 


which proves that the factorization theorem is violated at all order in coupling constant 
when the non-light-like Wilson line is used in QCD in the general Coulomb gauge where v^ 
is the non-light-like four velocity. 

Hence we have proved in this paper, by using path integral method of QCD, that the soft- 
collinear divergences are factorized at all order in coupling constant when light-like Wilson 
line is used in the definition of the parton distribution function and fragmentation function 
at high energy colliders but the soft-collinear divergences are not factorized when the non¬ 
light-like Wilson line is used in the definition of the (transverse momentum dependent) 
parton distribution function and fragmentation function at high energy colliders. This is 
in conformation with the finding in [7), 8] which used the diagrammatic method of QCD 
to prove factorization theorem in NRQCD heavy quarkonium production in case of light¬ 
like Wilson line [7] and violation of factorization theorem in NRQCD heavy quarkonium 
production in case of non-light-like Wilson line js]. In case of massive Wilson line in QCD 
the color transfer occurs and the factorization breaks down |8|. 

nr 

Hence we find that the factorization theorem is violated in all the previous studies [l|, 6] 
which used the non-light-like Wilson line in the definition of the (transverse momentum 
dependent) parton distribution function and fragmentation function at high energy colliders. 


IX. CONCLUSIONS 

By using path integral formulation of QCD and QED we have proven that the factor¬ 
ization theorem is valid for light-like Wilson line but is not valid for non-light-like Wilson 
line. This conclusion is shown to be consistent with Ward identity and Grammer-Yennie 
approximation. Hence we have found that the factorization theorem is violated in all the 
previous studies which used the non-light-like Wilson line in the definition of the (trans¬ 
verse momentum dependent) parton distribution function and fragmentation function at 
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high energy colliders. 
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